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Abstract 

Consistent perturbation theory for thermodynamical quantities in type II 
superconductors in magnetic field at low temperatures is developed. It is 
complementary to the existing expansion valid at high temperatures. Magne- 
tization and specific heat are calculated to two loop order and compare well 
to existing Monte Carlo simulations and experiments. 
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Thermal fluctuations play much larger role in high Tc superconductors then in the low 
temperature ones because the Ginzburg parameter Gi characterizing fluctuations is much 
larger. In the presence of magnetic fleld the importance of fluctuations in high Tc supercon- 
ductors is further enhanced. Strong magnetic fleld effectively suppresses long wavelength 
fluctuations in direction perpendicular to the fleld reducing dimensionality of the fluctuations 
by two [|l|. Under these circumstances fluctuations influence various physical properties and 
even lead to new observable qualitative phenomena like vortex lattice melting into vortex 
liquid far below the mean fleld phase transition line. It is quite straightforward to system- 
atically account for the fluctuations effect on magnetization, speciflc heat or conductivity 
perturbatively above the mean fleld transition line using Ginzburg - Landau description . 
However it proved to be extremely difficult to develop a quantitative theory in the inter- 
esting region below this line, even neglecting fluctuation of magnetic fleld and within the 
lowest Landau level (LLL) approximation. 

To approach the region below the mean fleld transition line T < Tmf{H) Thouless 
proposed a perturbative approach around homogeneous (liquid) state was in which all the 
"bubble" diagrams are resummed. The series provide accurate results at high temperatures, 
but for LLL dimensionless temperature = (T — Tmf{H))/{THy/^ < —2 become inap- 
plicable (see dotted lines on Fig. 2, 3 which represent successive approximants). Generally 
attempts to extend the theory to lower temperature by Fade extrapolation were not success- 
ful and require additional external information about the low temperatures . Alternative, 
more direct approach to low temperature fluctuations physics might have been to start from 
the Abrokosov solution at zero temperature and then take into account perturbatively devia- 
tions from this inhomogeneous solution. Experimentally it is reasonable since, for example, 
speciflc heat at low temperatures is a smooth function and the fluctuations contribution 
experimentally is quite small. This contrasts sharply with theoretical expectations. 

Long time ago Eilenberger calculated spectrum of harmonic excitations of the triangular 
vortex lattice Subsequently Maki and Takayama noted that the gapless mode is softer 
then the usual Goldstone mode expected as a result of spontaneous breaking of translational 
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invariance. The propagator for the "phase" excitations behaves as l/(fc^ + c(/c^ + ky)). The 
influence of this unexpected additional " softness" apparently goes even beyond enhancement 
of the contribution of fluctuations at leading order. It leads to disastrous infrared divergen- 
cies at higher orders rendering the perturbation theory around the vortex state doubtful. For 
example contributions to energy depicted on Fig. lA and ID are respectively log^(L) and 
divergent {L being an IR cutoff) and at higher orders divergencies get worse. Also quali- 
tatively one argues (in a way similar to that used frequently to understand the Mermin - 
Wagner theorem |^) that lower critical dimensionality for melting of the Abrikosov lattice is 
D = 3 and consequently vortex lattice in clean materials exists in the thermodynamic limit 
only at T = 0. One therefore tends to think that nonperturbative effects are so important 
that such a perturbation theory should be abandoned and it was abandoned. However 
a closer look at the diagrams like Fig. ID (see some details below) reveals that in fact one 
encounters actually only logarithmic divergencies. This makes the divergencies similar to so 
called "spurious" divergencies in the theory of critical phenomena with broken continuous 



symmetry. In that case one can prove |]IU| that they exactly cancel at each order provided 
we are calculating a symmetric quantity. 

In this note I show that all the IR divergencies in free energy or other quantities invariant 
under translations cancel to the two loop order. I calculate magnetization and speciflc heat 
to this order, interpolate with existing high temperature expansion and compare with Monte 
Carlo (MC) simulation [jlT| of the same system and experiments. Qualitatively physics of 
fluctuating D = 3 GL model in magnetic fleld turns out to be similar to that of spin systems 
(or scalar flelds) in D = 2 possessing a continuous symmetry. In particular, although within 
perturbation theory in thermodynamic limit the ordered phase (solid) exists only at T = 0, 
at low temperatures liquid differs very little in most aspects from solid. One can effectively 
use properly modifled perturbation theory to quantitatively study various properties of the 
vortex liquid phase. 

The GL free energy is 
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Here A = {—By, 0) describes a nonfluctuating almost constant magnetic field in c direc- 
tion. Within the LLL approximation ip can be expanded in a basis of quasimomentum k 
eigenfuntions 
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Unit of length will be the magnetic length Ifj = ^Jch/eH and ca = y 47r/v^ is the lattice 
spacing. The k = component ^o{x) = V'(a;) is "a vacuum" with its VEV denoted by v. 
The integration is over Brillouin zone. Instead of one complex field two real fields O and 
A were introduced. They are somewhat analogous to acoustic and optical phonons in usual 
solids with some peculiarities due to strong magnetic field studied in detail by Moore 0. For 
example the A mode corresponds to shear of the two dimensional lattice. Substituting eq.(|D 
into free energy, quadratic terms in fields define propagators, while cubic and quartic are 
interactions. The phase factors containing a function 7k = ^^ip*{x)ip*{x)({)\i{x)ip-\^{x) are 
introduced in order to diagonalize the resulting quadratic part PQ^{k)OlOk + PX^{k)AlAk, 
where PQ\{li,kz) = 2a + 26f^(2/3k ± |7k|) + k"^ (to simplify intermediate expressions an 
isotropic case rriab = rric is considered, results are generalized later). Functions 7k = A(— k, k) 
and /3k = '^*{x)v{x)vti^)Vi<.{^) = -^(O; k) as well as all the three and four leg vertices can 
be expressed via single function of two quasimomenta 
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. If the fluctuations were absent 



where (3^, ^ (^|^)' ^^^^ /«^*(_)^- exp [-^(/^ + m)^ 
the expectation value = would minimize Gq = —av"^ + \(3av^ where [3a = Pqq = 1-16. 
The propagators entering Feynman diagrams therefore are: 



Expanding around k = using explicit expressions for 7k and /3k one observes that constant 



and the k"^ terms vanish, so that the (only) leading quartic term is M^(k) = ^^|k|^. 



At one loop level the fluctuation contribution to the free energy is: 

Gi = l^^, 1^ ^ {log[Pom + \og[PA{k)]} (6) 
One should minimize Go + Gi with respect to v leading to the correction to its value: 
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Due to additional softness of the A mode the second "bubble" integral diverges loga- 
rithmically in the infrared. This means that for the infinite cutoff fluctuations destroy the 
inhomogeneous ground state, namely the state with lowest energy is a homogeneous liquid 
|r5|. Since the divergence is logarithmic we are at lower critical dimensionality in which 



an analog of Mermin - Wagner theorem [§] is applicable. It however does not necessarily 
means that perturbation theory starting from ordered ground state is inapplicable. The 
way to proceed in these situations have been found while considering simpler models like 
y?^ model F = |(V<y9j)^ + V{ipf) in D = 2 with number of components larger then 1, say 



i = 1,2 |T2|. Considering statistical sum, one first integrates exactly zero modes existing 
due to continuous symmetry (translations in our case) and then develops a perturbation 
theory via steepest descent method for the rest of the variables. When the zero mode is 
taken out, there appears a single configuration with lowest energy and steepest descent is 
well defined. For invariant quantities like energy this procedure simplifies: one actually can 
forget for a moment about integration over zero mode and proceed with the calculation as if 
it is done in the ordered phase. The invariance of the quantities ensures that the zero mode 



integration trivially factorizes. This is no longer true for noninvariant quantities for which 



the machinery of "collective coordinates method" should be used ||14|| . 

To the two loop order one gets several classes of diagrams, see Fig.l. The leading order 
propagators are denoted by dashed and solid lines for the "supersoft" A and "hard" B modes 
respectively. The naively most divergent diagram Fig. ID actually converges. To see this 
one writes explicitly its expression in terms of function A 
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The integrals over and Qz can be explicitly performed using a formula 
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/p /q -^"(q^ P) p2q2|q+p|^ p2+q2||q+p|2 , IS determined by the asymptotics of /^(q.p) as both 
p and q approach zero. If //) ~ 1, it would diverge as L^. However the vertex is "supersoft" 
at small quasimomenta ~ according to eq.(0), so that expansion of /£i(q, p) starts from 
terms quartic in p and q and there is no singularity at the origin. This goes beyond the 
usual "softness" of interactions of the Goldstone modes (~ p). Nonleading divergences can 
be found by analyzing contributions coming from three regions on which one of the line 
momenta p, q or p + q vanishes. The corresponding expressions are ~ -^d(I) k^MAO)'^ ' 
with l}j = 0, /|, = — /?i|7i| and /|) = —^^ + /3i|7i| respectively. Here k denotes an IR 
divergent momentum while integration over 1 is nonsingular. Although the second and the 
third contributions are divergent their sum is convergent. 

Standard methods similar to one used above can be applied to evaluate IR divergencies 
of other superficially less divergent diagrams. There are no power divergencies - only log^ L 



and logL.The results are ^ [ ^ [ + , ^ f ^ [ '^'^l^'f^'' 

'=' V27r Jp Ma{p) Jq [ Ma{ci) A/o(q)J ' V227r Jp Ma(p) Jq MA(q) 

7^ Ip Jljjp) /q '^Molnt ^'^^ diagrams Figl a, b and e respectively. In addition to direct 
contributions from G2 Fig.l there is also a "correction term" due to correction in the value 
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of V from eq.(|^) inserted into the lower order contributions to free energy Gq and Gi. It's 
divergent part is /-^ [ ult^)^^ - ^I&if^J • 2°*^ the leading divergencies 

log^L and the next to leading ones logL cancel between the four contributions. Similar 
cancellations of all the logarithmic IR divergencies occur in scalar models with Goldstone 
bosons in D = 2 and D = 3 (where the divergencies are known as "spurious"). 

The finite result for the Gibbs free energy to two loops (finite parts of the integrals were 
calculated numerically) is restoring the original units and reintroducing the asymmetry 
rric ^ mah- 
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where numerical values of the coefficients are Ci = 3.16 and C2 = 7.5. Dimensionless entropy 
(LLL scaled magnetization) 

and specific heat normalized to the mean field value 
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for successive partial sums are plotted on Fig. 2 and 3 (dashed lines). Qualitative they 
are in accord with numerous experiments and MC simulations |Tl|. At low temperature 
magnetization is a bit larger then the mean field's one, while dimensionless specific heat 
characteristically grows before dropping fast around qt = —5. To make more detailed 
comparison, I interpolated between results of low temperature expansion and these of high 
temperature expansion using the following rational form for free energy in terms of often 
used variable x defined imphcitly by x = y^, ax = 4(2?/)^/'^ (1 — 1/8?/^): 

.,3iW_±^^ (12) 

The coefficients were constrained from both low and high temperature sides. It has been 
already noted [§] that constraining from both sides the Fade approximants just by the 



first term at low energy improves otherwise unsatisfactory magnetization and specific heat. 

Adding two more terms on the low temperature end makes it very close to the MC results 
(stars, triangles and diamonds correspond to IT, 2T and 5T results for YBCO). I used just 
three leading terms in high temperature expansion shown on Fig. 2 and 3 by dotted lines. 
Using more terms does not modify significantly the result. Although magnetization curve 
eq.(|l2D agrees with that of ref. [|l^, the specific heat is not. 

To summarize, it is established up to order of two loops that perturbation theory around 
Abrikosov lattice is consistent. All the IR divergencies cancel due to soft interactions of the 
soft mode. Perturbative results as well as interpolation with the high temperature expansion 
agree very well with the direct MC simulation. 

Now I comment on range of validity of the perturbative results and nonperturbative 
effects. As can be seen from Fig. 2 and 3 the range of validity of the low temperature 
expansion presented in this paper is below = —10, while that of the high temperature 
expansion is above = —2. Both exclude the range in which small magnetization jump 
(not seen on Fig.2's scale) due to vortex melting is seen experimentally and in the numerical 
simulation. Since the MC simulation is the only systematic tool available in the intermediate 



region (the theory of Tesanovic et al ||T5|] captures major (98%) contribution, but does dot 



treat the small (2%) effect including melting), one might have two possibilities to discuss 
such a singularity within the present framework. One possibility is that the jump is due 
to finite size effects and disappears in the thermodynamic limit (value of the cutoff in the 
simulation is only L ~ 25). Another is that some nonperturbative effects can stabilize the 
vortex lattice. Quantitative comparison with experiments on YBCO was attempted in ref. 



rT|,|T5|. The present simple interpolation formula eq.([T^) works equally well. 
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Figure Captions 
Fig.l. 

Contributions to the free energy at the two loop leveL 
Fig. 2. 

Scaled magnetization defined in eq.([TD|). Dashed (dotted) lines are successive low (high) 
temperature approximants, while the solid line is the interpolation. Points are the MC 
results. 

Fig. 3 

Scaled specific heat defined in eg . (pA]) . Same notations as in Fig. 2. 
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